(N 
(N 






c^ 



MPI-PhT/2000-37 



Maximum lepton asymmetry from active — sterile 
neutrino oscillations in the Early Universe. 

R. Buras f\ 

'(^ • Max-Planck-Institut fiir Physik (Werner- Heisenherg-Institut) 

o: 

O ■ Fohringer Ring 6, 80805 Miinchen, Germany 

(N 

f^ , u. \ . ot:iiiiK.uZi rj 

CD ■ 
^yy I Max-Planck-Institut fiir Physik (Werner-Heisenherg-Institut) 



Fohringer Ring 6, 80805 Miinchen, Germany 
and 



\Q \ Institute of Nuclear Research of the Russian Academy of Sciences 

\o ' 

C^ • 60th October Anniversary Prospect 7a, Moscow 117312, Russia 

a^ 
o 
o 

O ' Abstract 



A large lepton asymmetry could be generated in the Early Universe by 
oscillations of active to sterile neutrinos with a small mixing angle sin 26 < 10~^. 



Q I The final order of magnitude of the lepton asymmetry rj is mainly determined 

1-^ ' by its growth in the last stage of evolution when the MSW resonance dominates 

the kinetic equations. In this paper we present a simple way of calculating the 
K^ \ maximum possible lepton asymmetry which can be created. Our results are 

;_( ■ in good agreement to previous calculations. Furthermore, we find that the 



growth of asymmetry does not obey any particular power law. We find that 
the maximum possible asymmetry at the freeze-out of the n/p ratio at T '^ 
IMeV strongly depends on the mass-squared difference 5m?: the asymmetry 
is negligible for \6m?\ <C 1 eV and reaches asymptotically large values for 
|5m2|^50eV2. 
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1 Introduction 

Neutrino oscillations between muon and tau neutrinos, which were found recently in 
the Super-Kamiokande experiment [IH , prove the existence of neutrino masses. This 
finding does not spoil the Standard Big Bang Nucleosynthesis (BBN) picture because 
experimentally interesting neutrino masses for all three flavors Ui, i = e,fi,T, are 
small rrii < 1 eV, in comparison to typical BBN temperatures T ~ 1 MeV. As for 
oscillations between active flavors, these do not play any role in BBN because all 
active flavors are equally populated. 

However, if sterile neutrinos exist, oscillations between active and sterile states 
could have significant consequences for BBN. One possibility is that the total energy 
density of the Universe is changed due to excited sterile neutrinos; together with a 
recent analysis of observational data [Q], which claims that the effective number of 
neutrinos can not exceed 3.2, this effect can be used to exclude some region of the 
{6m'^, sin 26) parameter plane e.g. for heavy sterile neutrinos 0. 

Another possible influence of sterile neutrinos on BBN arises from the fact that 
z/g-t's oscillations can create a large asymmetry between z/g and z/g which would have a 
direct impact on the n-p reactions and thus change the light element abundances. This 
case is interesting for small negative Sm^ and small sin 2^ < 10^^, and the important 
result will be the value of the electron neutrino asymmetry at the freeze-out of the 
n/p ratio at a temperature around 1 MeV. Actually, the magnitude of u^ asymmetry 
at T = 1 MeV is just an indicator for possible consequences on BBN. For accurate 
calculations one would need to use exact non-equilibrium distribution functions of z/g 
and z/g in the neutron-proton reactions. A general discussion of this subject can be 
found in Section 2.1 of [Q, and details in case of large asymmetry in 0. 

In early works with simplified calculations it was found that the asymmetry could 
not reach large final values @, |^. However, later it was found that the asymmetry 
can increase significantly and reach values (9(0.1) p|. This finding was confirmed 
by elaborate numerical calculations ^, |10], [ll]]. A main feature of these results was 



that the asymmetry grows according to the power law T ^. A contradictory result 
that the power law was T~^ was suggested in |]l2i, but was found to be erroneous in 
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In the present paper we investigate the evolution of the lepton asymmetry under 
the assumption that all active neutrinos passing the resonance are converted to sterile 
neutrinos, thus forming a Fermi-Dirac distribution in the sterile sector. This assump- 
tion is correct in the case of pure Mikheyev-Smirnov-Wolfenstein (MSW) transition 



151 , i.e. when one can neglect the collision-integral terms. Although the assumption 
of full conversion from active to sterile states may be too strong for certain models, 
our result does present the maximum possible asymmetry. The advantage of our 
approach is that we do not solve any differential equations. Thus, our approach is 
extremely simple and safe. Furthermore, we consider the effects of various simplifi- 
cations and compare our results with previous works, in particular with 0. Finally, 
we investigate how the maximum amplitude of the asymmetry at the temperature 
1 MeV depends on initial conditions and 6m^, and compare it with a more detailed 
calculation, [|T^ . 

2 Maximal asymmetry at given temperature 

We are considering neutrino oscillations in the Early Universe at temperatures around 
1 MeV < T < 100 MeV. In this epoch, the Universe is a homogeneous and isotropic 
plasma, and the expansion is described by the Hubble parameter. 
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where nipi = 1.22 x 10 MeV is the Planck mass. The total energy density is 

Ptot = ^9*T^ ' (2) 

where g^, = 10.75 is the number of effective degrees of freedom during the epoch of 
interest. We also suppose that T/T = —H, which is a good approximation for the 
models we are considering here. 



The neutrino oscillations are described by the density matrix formalism, see for 



instance JT^ • However, we will mainly be interested in the resonance condition, which 
for Us, ^->- Vs oscillations is given by p 



'-^^ --v^-~ [^c^no.^ + c-^) . (3) 

where a = e, fi,T specifies the type of active neutrino, Sm"^ = ml — m^ is the mass 
squared difference, 9 is the mixing angle, and V^^ is the effective potential. The 
"=F" signs refer to neutrinos and anti-neutrinos, respectively. Furthermore, Gp = 
1.166 X 10~^ GeV^^ is the Fermi coupling constant, a = 1/137 is the fine structure 
constant, Ci = 0.345, C2 = 0.61 for Ue-i^s mixing, 6*2'^ = 0.17 for v^^t-i-'s and T is the 
temperature of the plasma. Sometimes we will use the notation (5m^y = 5m^ /eV"^. 
Finally, the effective asymmetry is r^ = 2ri„^ + riQ, where 7],^^ is the neutrino asymmetry 
of the active species a, and r^o is some small fermion asymmetry from all other fermion 
species. For our analysis it is only interesting that rjQ is of order of the baryon 
asymmetry, 10~^°. 

Because the temperature is much larger then the neutrino masses, not only neu- 
trino momentum and energy are equal, p = E, but also that a certain neutrino mode 
evolves with constant y = p/T. We will therefore use this dimensionless variable 
instead of momentum and energy, so p = E = yT. 

We briefly discuss the evolution of the neutrino densities for small mixing angles, 
sin 26' < 10~^, looking at u^ for definiteness. Let us first sketch the reason why the 
lepton asymmetry can change at all. Due to the mixing of u^ with z/g, the number 
densities of the z/g and z/g can alter. The sterile neutrinos that are thus created have 
no effect on rj. Now if r^ 7^ initially, the resonance condition for the z/g and z^e 
will be different. Therefore their number densities will change differently, thereby 
changing rj. Thus neutrino mixing works as a back-reaction on rj. The actual change 
in asymmetry and density takes place close to the resonance momentum, which is 
defined by Eq. (|), 

yres = j[h^y + k,^±h^, (4) 



^1 = o^e^ /. and k2 = ^^,^2 /. • (5) 



where 

_C7i _ |^m^|cos2g 

There are three different stages of lepton asymmetry evolution. During the first 
stage the asymmetry is so small that the ki terms can be neglected in Eq. (^). The 
neutrinos and the anti-neutrinos simultaneously pass through the resonance. Fur- 
thermore, any initial asymmetry decreases due to the negative back-reaction. 

The second stage is the stage of exponential growth of the lepton asymmetry, which 
occurs when the back-reaction becomes positive. Here, the resonance conditions for 
z/ and u are rapidly driven apart. For definiteness, we will assume that the created 
asymmetry is positive, r/ > 0. Then ^rcs ioT the u increases to very high momenta 
where the Fermi distribution is negligible, the resonance momentum for z/ decreases 
to small momenta ^rcs < 1. At the same time, the asymmetry grows many orders of 
magnitude and reaches values of order of 10~^. 

In the third stage the lepton asymmetry grows at a slower rate. There has been 
a lot of controversy on this subject, and it has often be proposed that the lepton 
asymmetry would grow according to an approximate power law, rj oc T^", where 
a = 1-4. 

We will now show that any power law is only an approximation to the true asym- 
metry growth. Our considerations start at the end of the exponential regime, which 
happens at temperature Tj. The important variables are the asymmetry i]i{Ti) and 
the resonance momentum yi{Ti,rii) at this time. For definiteness, we will assume 
rji > 0, so that we can safely neglect z/ oscillations (with r/j < we could neglect u). 

The resonance momentum is governed by two effects: it decreases for increasing 
rj and increases for decreasing T. At the end of the exponential regime, these two 
effects compensate each other, so that at first the resonance momentum changes very 
slowly. This plus the fact that collisional damping is small may let us assume that 
the momentum modes passing the resonance in the regime of slow asymmetry growth 
will fully excite the sterile sector. For negligible collisional damping, this is exactly 



the MSW effect [^. The full excitation of the sterile sector at resonance has also 
been stated in 0. 

We can now calculate the asymmetry for a given temperature T < Ti'. 

viT) = Vi + ^T7^ J dyy^feqiy, P) ■ (6) 

Here, /eq is the equilibrium Fermi distribution function, given by 

/eq(l/,/2) = ^^_^^^ , (7) 

where jl = fi/T is the dimensionless chemical potential for the active neutrinos. For 
1] <^ 1 and large temperatures T ^ 1 MeV, the chemical potential is approximately 
/i ^ — l.Sr/. Therefore, we can safely neglect it. We will discuss its influence for large 
rj later on. We should mention that rji has its origin in an already weakly excited 
sterile sector at momenta < y < 2. However, the correction to Eq. @ from this 
excitation is negligible. 

We can now insert Eq. (^) into Eq. (|^) and solve the equation for rj{T). The 
resulting equation only depends on the input parameters rji and yt, and on the neu- 
trino parameter Sm^. Interestingly, the equation is independent of the mixing angle 
sin 26', a fact which stems from the full transition assumption. The result has another 
striking feature: even if the sterile sector is not fully excited during resonance, which 
happens if the resonance momentum changes too fast and/or the coUisional damping 
is important, our result still gives the maximally possible asymmetry at temperature 
T. 

For small momenta, the Fermi-Dirac distribution can be approximated by /eq = 
1/2. Then Eq. @ simplifies to 

"^(^^ - "^^ = Y^) {y'e^^'^^ - y') ■ (8) 

Furthermore, for large asymmetries and small temperatures, the ^2 term in Eq. (1) 
is much smaller than the ki term, so that 

k2 



2ki7]T^ 
6 



(9) 



Then 

^' + ^' (12^^^' - ^^) = 12^ (2^) • ^'°) 

This equation can be easily solved in two limits: when rj <^ ( i2C(3) y^ ~ Vi) ^^ can 
neglect the first term on the LHS so that we obtain a. rj cc T~^ power law. However, 
when r] ^ ( i2C(3) ^»^ ~ ^^^ ^^^ solution is rj oc T~^. Thus the power law will slowly 
change in time. 

This is also what we expect from the following considerations: the asymmetry 
is created close to the resonance momentum. However, when the sterile sector is 
fully excited at this momentum, the asymmetry can no longer be created. Naturally, 
when 7] does not change, the resonance momentum will increase due to the decreasing 
temperature. Thus we expect a slowly increasing resonance momentum. From Eq. (^) 
we can then deduce that r] increases more slowly than a. rj cc T~^ power law. 

We can easily derive an equation similar to Eq. (14) of [^ by differentiating Eq. (^) 
in the temperature T: 

The only difference between our Eq. (^) and Eq. (14) of [Q is in the usually negligible 
term proportional to f^^. Furthermore, Eq. (18) of 0, which can be derived from 
their Eq. (14) directly, is in our notation 



drj _ ~^2a^f'=^ 



(12) 



where we have used the simplification Eq. @. 

For illustration, we have plotted the evolution of the "power" a which can be 
derived by the formula 

« = -^ («) 

rj al 
in Fig. |l]. Here, we have already included the effect from the exact Fermi-Dirac 

distribution, which becomes important at low temperatures. We have plotted the 

evolution for 5m'^ = 1 eV^ and for different values of yi. Note that typical values are 
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Figure 1: Power-law index a as a function of temperature for various values of the ini- 
tial resonance momentum ?/j. We have used the simplified condition for the resonance 
momentum and have neglected the initial asymmetry rji. Note that e,g. r^j ~ 6 x 10~^ 
would substantially change the line for t/j < 0.1. 
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Figure 2: Power-law index a as a function of temperature for various values of the 
initial resonance momentum yi. Here we have used the exact resonance momentum 
condition. The solid lines represent the case rji = 0, whereas the dashed lines show 
the effect from non-zero rji, in this case rji = Q x 10~®. The dots show starting-points, 
i.e. Ti and the corresponding power aj (for r/j = 0, Tj = oo). 



yi = 0.1-0.2. Furthermore, we have set rji = 0. Typical values would be rji = 10^^- 
10^^, but as we can see from Eq. ([T0|), such values are negligible in comparison 
with the yf term. Anyhow, we can simulate the effect of rji by changing yi, see 
Eq. (10). Interestingly, the effect from 6m'^ can be compensated exactly by rescaling 



the temperature: 

Tsm^ = {6m'Y/' T,^,^,^y2. (14) 

To be still more precise, we need to use the exact Eq. (Q) instead of Eq. (|^). We 
have plotted the results in Fig. || in analogy to Fig. |l[ We see that both the approx- 
imate and the exact resonance condition give the same results at low temperatures, 
while at high temperatures the results are very different; the exact result shows a 
larger power-law behavior. This means that there is no real power-law solution in 
the MSW dominated region, at least under the assumption that transition is com- 
plete. We should mention that for v^^r-^s oscillations, the correction due to the exact 
resonance condition is somewhat weaker, since C2 is smaller by a factor 4. 

We should note that using Eq. (P) instead of Eq. (H) has no influence on the later 
evolution of the asymmetry, i.e. Eq. (^) is a good approximation. This is because the 
apparently larger increase of r] for given Tj, rji is compensated by a smaller yi when 
using Eq. (|). 

In Fig. 1^ we also show the effect of a non-zero initial asymmetry. This effect is 
important only for small initial values of the resonance momentum |/res < 0.1. 

3 Asymmetry at the freeze-out of the n/p ratio. 

If we want to investigate the influence of the neutrino asymmetry on BBN, we need 
to check its value at the temperature T = 1 MeV when the neutron-proton reactions 
decouple and the n/p ratio freezes out. In Fig. ^ we present the dependence of the 
maximum possible asymmetry at the temperature 1 MeV from |(5mgY|. The solid line 
shows the maximum asymmetry for the initial conditions yi = and rji = 0. Our 
result is in good agreement with previous extensive calculations |TB]. Naturally, we 
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Figure 3: The maximal possible asymmetry t] at the freeze-out temperature of the 
n-p interactions, T = 1 MeV, as a function of |(5m^v|, for rji = i/i = (solid line). The 
initial asymmetry rji increases the final asymmetry for small If^mgyl (short dashed 
line for rji = 10^"^); the momentum i/i decreases the asymmetry for small l^mgyl (long 
dashed hue is for yi = 0.2). 
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Figure 4: The resonance momentum |/res at the freeze-out temperature of the n-p 
interactions, T = 1 MeV, as a function of l^m^yl. The initial asymmetry rji and 
momentum yi only have significant effects for very low l^m^yl; we have neglected 
them in this plot. 
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do not expect any dependence on sin 29, since we have assumed full MSW transition. 
We should therefore stress that our result is only a good approximation in regions of 
the parameter space (5m^,sin26') where there is an effective exponential increase of 
the asymmetry, and if there is full MSW transition afterwards. 

From Fig. ^we can easily extract that for l^m^yl <^ 1 the asymmetry will remain 
small, while for l^wigyl ^ 1 it can reach its asymptotic value 3/8, which was found in 
0] (the factor 2 difference with our 6/8 on Fig. |^ stems from the relation r] = 21]^ J. 
For l^^vl ^ -'-' ^^^ asymmetry at temperature 1 MeV can be fitted according to the 
simple relation 

ri = O.l X \5mlyf^ (15) 



in good agreement with the result of [0, rj = 0.1 x \6m^ 



2 12/3 



We will shortly discuss minor effects on the result. The short dashed line in Fig. ^ 
shows the effect of a non-zero initial asymmetry (for rji = 10~^ in this particular case). 
A non-zero initial asymmetry will increase the asymmetry at T = 1 MeV for small 
|5?7i^|. A non-zero momentum yi will work in the opposite direction and decrease the 
final asymmetry for small I^T^igyl (long dashed line is for t/i = 0.2). Note that for very 
small |5mgvl ^ 10^^ the exponential regime is to be expected at T ~ 1 MeV so that 
our results become invalid. Finally, the chemical potential plays an important role 
when the asymmetry is large. However, its effect is to exponentially damp the number 
density of the z/g (in the case ?7 > 0), thus decreasing the asymmetry. Therefore, our 
result remains the maximum possible asymmetry. We have estimated the decrease of 
the asymmetry due to the chemical potential to be less than a factor of 2. 

For illustration. Fig. ^ shows the dependence of the dimensionless resonance mo- 
mentum i/res = Pies/T at the temperature T = 1 MeV from the parameter I^TTigyl. 
We see that at If^TTigyl <^ 1 the resonance momentum is small, while for |5?7igy| > 10 
it is already on the tail of the distribution function. This explains the values of the 
asymmetry presented in Fig. 13. 
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4 Conclusions 

In this work we have found the maximum lepton asymmetry which can be created at 
a given temperature due to active-sterile neutrino oscillations in the Early Universe. 
We have assumed full transitions of active to sterile neutrinos due to MSW resonance 
to derive a very simple relation between the lepton asymmetry and the temperature. 
From this relation, we found that the evolution of the lepton asymmetry does not 
evolve according to a specific power law. 

We have also investigated on the maximum asymmetry at the temperature 1 MeV 
when the n/p ratio freezes out. We found that the value of the asymmetry at this 
time depends on 5m?, but not on sin 2^ in concordance with |]T6|. Furthermore, the 
initial settings described by the initial resonance momentum yi and asymmetry rji 
have only an influence for very small 5m?. This means that details of the previous 
evolution are not very important for the value of the asymmetry at T = 1 MeV. In 
the parameter region which is related to a large final asymmetry, 5m1y > 10, our 
asymmetry values are in good agreement with those of |^. 

We stress again that our results depend on the assumption of full transitions of 
active to sterile neutrinos. If in some part of the parameter space the transition is 
not total, then the asymmetry can be somewhat smaller. Also, for large asymmetries, 
the exact value of rj may be smaller due to effects from the chemical potential. 

Acknowledgment 

We are grateful to G. Raffelt for helpful comments and reading the manuscript. We 
also would like to thank S. Hansen, R. Foot and R.R. Volkas for useful comments. 

This work was partly supported by the Deutsche Forschungsgemeinschaft under 
grant No. SFB 375 and in part by INTAS grant 1 A- 1065. 



12 



References 

[1] Super-Kamiokande Collaboration, Y. Fukuda et ai, Phys. Rev. Lett. 81 (1998) 
1562; ibid. 82 (1999) 2644; ibid. 82 (1999) 5194; Phys. Lett. B467 (1999) 185. 

[2] D. Tytler, J. M. O'Meara, N. Suzuki and D. Lubin, to appear in Physica Scripta, 
|astro-ph/0001318 . 



[3] A.D. Dolgov, S.H. Hansen, G. Raffelt and D.V. Semikoz, |hep-pli/0008T^. 



[4] A.D. Dolgov, S.H. Hansen, G. Raffelt and D.V. Semikoz, Nucl. Phys. B580 
(2000) 331. 

[5] R. Foot, Phys. Rev. D61 (2000) 023516. 

[6] K. Enqvist, K. Kainulainen and J. Maalampi, Phys. Lett. B244 (1990) 186. 

[7] R. Barbieri and A. Dolgov, Nucl. Phys. B237 (1991) 742. 

R. Foot, M. Thomson and R.R. Volkas, Phys. Rev. D53 (1996) 5349. 

[9] R. Foot and R.R. Volkas, Phys. Rev. D56 (1997) 6653; Erratum ibid. D59 (1999) 
029901. 

[10] R. Foot, Astropart. Phys. 10 (1999) 253. 

[11] R Di Bari and R. Foot, Phys. Rev. D61 (2000) 105012. 

[12] A.D. Dolgov, S.H. Hansen, S. Pastor and D.V. Semikoz, Astropart. Phys. 14 
(2000) 79. 



[13] R. Buras, D.V. Semikoz, |hep-ph/0008263 



[14] P. Di Bari, R. Foot, R. R. Volkas, Y. Y. Y. Wong, |hep-ph/00082l5 



[15] L. Wolfenstein, Phys. Rev. D17 (1978) 2369; ibid. D17 (1979) 2634; 
S.P. Mikheyev and A. Yu. Smirnov, Nuov. Cim. C9 (1986) 17. 

13 



[16] X. Shi, G.M. Fuller, K. Abazajian, Phys. Rev. D60 (1999) 063002. 

[17] A.D. Dolgov, Sov. J. Nucl. Phys. 33 (1981) 700. 

[18] D. Notzold and G. Raffelt, Nucl. Phys. B307 (1988) 924. 



14 



